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Abstract 

We present the full implementation of a room-temperature gravitational wave bar detector 
equipped with an opto-mechanical readout. The mechanical vibrations are read by a Fabry-Perot 
interferometer whose length changes are compared with a stable reference optical cavity by means 
of a resonant laser. The detector performance is completely characterized in terms of spectral sen- 
sitivity and statistical properties of the fluctuations in the system output signal. The new kind of 
readout technique allows for wide-band detection sensitivity and we can accurately test the model 
of the coupled oscillators for thermal noise. Our results are very promising in view of cryogenic 
operation and represent an important step towards significant improvements in the performance 
of massive gravitational wave detectors. 
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I. INTRODUCTION 



The direct observation of gravitational waves (GWs) is one of the most challenging tasks 
for experimental physics. The effort devoted to this goal started in the 1960 's, based on 
theoretical predictions of the expected signal, which are now considered as very optimistic [|1| . 
The experimental strategy for detecting GW signals was based initially on massive acoustic 
detectors and nowadays exploits also long baseline interferometers 0. In particular, the 
former are the most sensitive GW detectors presently in activity [Q and they offer interesting 
possibihties for future advanced versions ^ . 

Cryogenic bar acoustic detectors are currently equipped with capacitive or inductive 
transducers followed by SQUID amplifiers or by a microwave resonant cavity. The sensitivity 
is presently limited by the amplification stage that operates ~ 10^ times above the standard 
quantum limit 0, giving a bandwidth of few Hz around the two mechanical vibration modes 
of the coupled oscillators system. 

The possibility of using optical techniques for the readout of bar vibrations was early 
considered by Drever 0. Kulagin et al. p theoretically studied the possibilities of an 
optical readout system for a Weber bar with resonant mechanical transformer. This idea 



was developed by Richard [jTO|, who designed and theoretically investigated in detail such 
a system |TT[. Richard and coworkers also tested at room temperature an opto-mechanical 
transducer made by a Fabry-Perot cavity installed on a double oscillator, observing a rms 
displacement noise consistent with the calculated thermal fluctuations . 

In spite of this effort, GW bar detectors have never been equipped with optical readout 
and the real performance of such a system has never been experimentally verified. On 
the other hand, the technology has very advanced during the last years in the fields of 
laser stabilization and fabrication of optical components. As a consequence, the expected 
characteristics of an optical readout system are even more promising and could lead to a 
major breakthrough of GW massive detectors. 



In the framework of the AURIGA collaboration we are developing a complete optical 



readout system for ultra-cryogenic bar detectors [|14|. The basic idea is to form a high- 
finesse Fabry-Perot cavity between the bar and a resonant mechanical transducer and then 
to compare the length of this optical resonator, possibly carrying a GW signal, with that of 
a stable reference cavity by means of a resonant laser. 
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In this work we present the first room temperature GW bar detector operating with an 
optical readout system. The detector performance is completely characterized in terms of 
spectral sensitivity and statistical properties of the fluctuations in the system output signal. 
Our apparatus represents a wide-band (several tens of Hz), GW acoustic detector which is 
limited by thermal noise at least in the frequency range of highest sensitivity. Thanks to 
this property, we can accurately study the output spectrum of the thermal noise and we 
show that the description of the bar and transducer coupled oscillators cannot be given in 
terms of de-coupled normal modes, as usually assumed. 

The outline of this paper is the following. In Section II we describe the experimental 
apparatus, and in particular the readout system. In Sec. Ill A we present the mechanical 
characteristics and the displacement noise of the detector, as deduced from the output 



signal. The effect of thermal noise in the output spectrum is analyzed in Sec. |III B| . Then 
we investigate the statistical properties of the output fluctuations, performing the accurate 
analysis necessary to characterize a GW detector (Sec. Ill G ). Finally in Sec. Ill D| we 
estimate the sensitivity of the bar as GW detector. 



II. THE READOUT SYSTEM 



A schematic drawing of the room-temperature GW detector equipped with the opto- 
mechanical readout is shown in Fig. |l|. The mechanical vibrations of the bar are amplified 
by a coupled mechanical oscillator and transformed into length changes of an optical cavity, 
hereafter called transducer cavity (TG). The readout system is composed of a laser source, 
frequency stabilized to a reference cavity (RG), a set of optical fibers and components to 
convey the radiation to TG, the opto-electronics for signal detection and elaboration. 

The bar is a 3 m long, 2300 kg mass cylinder made of A15056. Its first longitudinal 
vibration mode, useful for GW detection, resonates at 875 Hz, when no load is applied. The 
measured resonance frequency varies as —0.38 Hz/kg with the amount of non-resonant mass 
applied to the bar end faces, giving a calculated frequency of z/b = 866 Hz when the bar 
is operated with the full readout system. The mechanical quality factor of the resonance 
is 1.8 X 10^, measured from the decay time of the excited vibration. The bar is kept in 
a vacuum chamber, placed a few meters apart from the optical table, and it is isolated 
from floor mechanical noise by a cascade of passive filters which achieve an overall vertical 
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FIG. 1: Experimental setup: the drawing is not to scale. 01: optical isolator; HWP: half-wave plate: 
QWP: quarter-wave plate; L: lens; EOMT^t; electro-optic modulator; P: polarizer; BS#: beam- 
splitter; PBS: polarizing beam-splitter; PD#: photodiode; RC: reference cavity; TC: transducer 
cavity. Dashed lines mark out active thermal stabilization. 

isolation of about —140 dB at the bar frequency, as sensed at the bar middle section. During 
the work here reported the vacuum system, composed by a roots pump backed by a rotary 
pump, was operated just for about 1 hour per day. 

The output beam of a commercial NdiYAG laser source, emitting 50 mW at 1.064 /im, 
passes through an optical isolator and two electro-optic modulators (EOMs) enclosed in a 
thermally stabilized box. The first EOM is used with a polarizer in an amplitude stabilization 
loop. The purpose of this noise-eater, described in detail in Ref. |T5|, is reducing the effect 



of the back-action on the transducer, but it is not relevant for the work here reported. 
The second EOM is a resonant modulator working at 13.3 MHz which accomplishes phase 
modulation with a depth of about 1 rad. A first beam-splitter (BSl) transmits 20% of the 
radiation for the noise eater, while the refiected beam is directed towards a second 70% 
transmission beam-splitter (BS2). The refiected beam, after an optical isolator, is coupled 
to a single-mode polarization-maintaining fiber and arrives to a 135 x 340 mm^ Al plate 
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anchored to the bar middle section. 

The optical fiber is formed by joining 4 patchcords with FC/PC connectors, for a total 
length of 13 m, and includes a homemade vacuum feedthrough. The two fiber ends have 
pigtailed collimators, with anti-refiection coating. The overall power transmission of the 
fiber assembly is about 50%. 

The 1.5 mW coUimated beam transmitted by the fiber passes through an optical circu- 
lator, formed by a polarizing beam-splitter and a quarter-wave plate, and a telescope to 
properly couple the radiation to the TEMqo mode of TC. Four tilting mirrors send the beam 
towards the transducer cavity on the bar end face. The beam refiected by TC, after the 
circulator, is detected by a photodiode (PDl). The TC is a 6 mm long Fabry-Perot cavity, 
with a finesse of 28000, formed by an input concave mirror (radius of curvature 1 m, diameter 
0.5") glued to a support fixed to the bar, and a fiat back mirror (diameter 0.5") glued to the 
oscillating mass of the mechanical transducer. This resonant transducer is machined from a 
single piece of A15056 and it is composed of a thin circular plate loaded by a central 1.25 kg 
inert mass. The resonant frequency of the first drum mode is about 882 Hz, according to the 
measurements described in Sec. ill A . This resonator was designed for a previous version of 
the readout system and it is described in Ref. 



The beam transmitted by BS2, after an optical isolator and mode- matching lenses, is sent 
to a 110 mm long Fabry-Perot reference cavity that has a finesse of 44000. RC is formed by 
an Invar spacer with a couple of mirrors similar to the ones of the TC. The input fiat mirror 
is glued on a piezoelectric actuator (PZT), which allows the tuning of the cavity length. 
The light power impinging on the cavity is about 4.5 mW. The cavity is kept in a vacuum 
chamber whose temperature is actively stabilized at about 34°C within 0.1°C. The beam 
refiected by this cavity is detected by a second photodiode (PD2) after an optical circulator. 

The power level impinging on PDl shows large variations: during the 42 hour period 
of continuous data acquisition, it varied by up to 80%, on a time-scale of typically a few 
hours. On the same period the power impinging on PD2 varied by less than 10%, with a 
longer timescale. The large variations sensed by PDl are due to polarization fiuctuations 
generated by drifts of the room temperature and originating from a non-perfect matching 
between the polarization axes of the fiber patchcords. They are turned into changes of the 
power impinging on TC by the optical circulator at the fiber output. 

The ac component of the signals coming from the two photodiodes is demodulated at 
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13.3 MHz and filtered, according to the Pound-Drever scheme |jT6[. The resulting signals 
are used as discriminator for frequency locking and analysis. The laser frequency is locked 
to a resonance peak of RC with a servo loop which has a unit gain frequency of 30 kHz. The 
loop gain is maximized in the frequency range between 600 Hz and 900 Hz, where it reaches 
120 dB, and allows to achieve an in-loop frequency noise level below the shot noise limit. 
A detailed description of the loop electronics and performance can be found in Ref. [IT? 



The resonant peak of RC is then superimposed to a resonance of TC by operating on the 
PZT actuator. The Pound-Drever signal from TC is used in a servo loop which drives the 
PZT through a low noise high-voltage amplifier. This low frequency servo loop has unity 
gain at about 1 Hz and allows the reference cavity to follow the resonance of TC in its 
thermal drifts. Due to the low loop bandwidth, the two cavities can be considered as free 
and independent in the frequency range of interest. 

The same error signal from TC is acquired and analyzed to extract information concerning 
the motion of the bar detector. 



III. EXPERIMENTAL RESULTS 

A. Mechanical system and noise spectrum 

The conversion of the detector output signal from voltage into length change of the 
transducer cavity is obtained using the slope of the corresponding error signal, measured 
with an accuracy of about 20%, and the known cavity length. The power spectral density 
Sxx of the displacement noise is shown in Fig. H, which corresponds to the average of one 
hour data. The peaks at 856 Hz and 892 Hz correspond to the frequencies v± of the two 
mechanical modes of the coupled oscillators system formed by the bar and the resonant 
transducer. For the coupled system the following relation holds: zz+z/^ = z/bZ/t, where Vi is 
the resonance frequency of the transducer. According to the value of z/b = 866 Hz calculated 
for the loaded bar we infer that the transducer resonance is at 882 Hz and thus detuned by 
+ 16 Hz with respect to the bar. A better coupling between the two oscillators is possible 
by adjusting the thickness of the transducer circular plate. 

We determined the mechanical quality factor Q of the '+' and ' — ' modes by measuring the 
decay time of a resonant sinusoidal excitation applied to the bar by means of a piezoelectric 



6 



N 




500 1000 
Frequency (Hz) 



T 1 1 1 r- 

1500 



FIG. 2: Power spectral density Sxx of the displacement noise. The two peaks at 856 Hz and 892 Hz 
corresponds to the two mechanical modes of coupled oscillators system formed by the bar with its 
first longitudinal mode and the transducer. Some other mechanical resonances are visible too. The 
sharp peaks at multiples of 50 Hz are due to the power line. 

actuator situated at the bar end face opposite the transducer. We got = 16600 and 
= 8700. The factor of 2 of difference in the Qs is explained by considering that, due to 
the frequency detuning between the two resonators, the minus mode is more influenced by 
the (high Q) bar while the plus mode by the (low Q) transducer. 

As it can be seen in Fig. |^, the peak spectral power in the modes exceeds by about 
45 dB the background noise. The wideband output noise of the optical readout system 
comes from the residual frequency fluctuations of the laser stabilized to RC and from the 
noise measured when the laser is far from the TC resonance. The laser frequency noise 
has been measured with respect to a stable Zerodur cavity, with an apparatus described in 
Ref. [18 1, and its effect gives a sensitivity limit as low as 2 x 10~^ Hz/\/Hz around 1 kHz. The 
far-from-resonance fluctuations are due to electronic noise, laser amplitude noise (including 
shot noise), and interference fringes. The overall effect gives a sensitivity limit of about 
0.2 Hz/v^Hz. The observed background, visible in Fig. 0, is about 20 dB higher than this 
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limit and exhibits a decreasing behavior versus frequency . 



We fitted the noise power spectral density of one hour output data with the function 
Ajv'^ between 75 Hz and 1775 Hz, neglecting only the interval around the '+' and '— ' modes 
and the resonance at 1519 Hz. We obtained e = 1.16 and A = A.2 x 10~^° m^Hz^Hz. The 
origin of this noise is unknown and is presently under investigation. We remark that a 1// 
frequency dependence in the noise spectrum is expected for the thermal noise of a mechanical 
oscillator with internal friction modeled as a constant imaginary part of the spring constant, 
in the frequency region below the oscillator resonance. 

B. Thermal noise of the coupled oscillators system 

A bar detector equipped with a resonant transducer is widely modeled as a system of 
two coupled harmonic oscillators and the output is often analyzed in terms of normal mode 
expansion. In particular, such assumption is used for studying the thermal noise of the 
system |^0[. On the other hand, it has been suggested and experimentally verified 



that this description may fail if inhomogeneously distributed losses occur. For our system 
the condition for the validity of the normal mode expansion can be written, according to 

as 



Ref. [21 



i^hQh = ^tQt ■ (1) 

In our case, the large difference between the quality factors of the modes '+' and '— ' can be 
brought back to a large difference in the effective Qs of the original oscillators Qb and Qt- 
Therefore, since the frequencies z/b and z/t are very similar, the condition of Eq. ([1|) is not 
satisfied. 

We call Xb, mt, and Xt, rrit the coordinate and effective mass of the bar and transducer 
oscillator respectively, and /b and /t the corresponding total driving force. Assuming that 
only viscous damping is present,^ the dynamics of the system is described by the equations 



^ We also considered the case of structural losses as dissipative mechanism, but no difference was found 
within the reported errors. 



8 



of motion 



Vb Vt 



+ 7=r (^t - a;b) + (^^t - a;b) = — 
Qt mt 



fh — ft 

mb 



where fj, = mt/m\^ and a;b,t = Svrz/b^f In the frequency domain the system can be written as 

(2) 



D 



(Fb - Ft)/mb 
Ft/mt 



with 



DM 



Qb 



"ft 



a; 



-^uo^ — I /i- 



Qt 



v 



— CJf — I 



9 r, . IX) UJt 

—ijj^ + u;^ + z 



where capital letters indicate Fourier transforms and i is the imaginary unit. 

To simplify our analysis, in the following we consider the transducer cavity length changes 
as determined exclusively by a motion of the transducer. This is justified as the amplitude 
of a bar displacement is amplified at the transducer by a factor equal to l/i//!, i.e., by 
a factor of about 30. We are thus interested to the noise power spectral density S^^xt of 
Xt which, if only stochastic thermal forces are present, can be written according to the 
Fluctuation-Dissipation Theorem ETI] as 

2k^T 



-Re{{iujT>)-\2} 



(3) 



where k-Q is the Boltzmann constant and T is the thermodynamic temperature. In order to 
account for the observed ~ 1// background, we add a phenomenological wide-band term: 

•S'wB 



(4) 



We fitted the output spectrum using the expression of Eq. in the frequency range 
between 780 Hz and 950 Hz. The fit allows to infer the ratios Qb/T", Qt/T and the resonant 
frequencies z/b, Ut of the uncoupled oscillators, the effective mass rrit of the transducer oscil- 
lator and the magnitude S'wB of the background noise. The effective mass mb = 1180 kg of 
the loaded bar resonator and the temperature T were kept constant during the fitting. We 



assumed that both oscillators are at the same thermodynamic temperature T of 296 K, as 
measured by a probe placed on the bar. This assumption seems reasonable at least for the 
bar as the inferred value for Qb agrees very well with the one measured independently for 
the unloaded bar (1.8 x 10^). The result of the fitting procedure is shown in Fig. ^ and the 
parameters are summarized in Tab. |. The two frequencies z/b and agree well with the 
values estimated in Sees. II and [111 A| . Also the scale factor S'wb of the background noise 
agrees with the value {2TTyA = 3.5 x 10~^^ m^Hz'^/Hz obtained in Sec. |III A| . As far as the 
transducer mass is concerned we notice that its effective mass deduced from the fit is greater 
than the mere 1.25 kg central mass. 

We also attempted to fit the same data with the prediction of the normal mode expansion, 
as shown in Fig. ^. It is evident that, while the experimental spectrum is in excellent agree- 
ment with the model of Eqs. and @, the normal mode expansion fails to describe our 
system: the presence of inhomogeneously distributed losses causes the random fiuctuations 
of the two initial oscillators to be correlated. The normal mode expansion overestimates the 
noise in between the two modes because it does not take into account such a correlation. 



C. Statistical analysis of the output fluctuations 

For a deep understanding of the detection system, it is important to investigate the 
statistical behavior of the noise before attempting any estimation of its magnitude. In fact, 
one needs to be sure that the noise under observation follows the laws predicted for the 

TABLE I: Results of the fit of the experimental spectrum with Eqs. (^) and (^), assuming T=296 K 
for both resonators. The quoted errors (two standard deviations) refer to the last significant digit. 



Parameter Fitted value units 

z/b 866.31 (3) Hz 

ut 882.30 (3) Hz 

Qb 1.8 X 10^ (4) 

Qt 6.60 X 10^ (4) 

SwB 3.5 X 10-29 (2) m^Hz-VHz 

mt 1.70 (2) kg 
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FIG. 3: a) Power spectral density of the displacement noise. Circles: experimental data; solid line: 
fit according to the two-oscillators model of Eqs. (|^ and dashed line: fit according to the 
normal mode expansion, b) The same as in a), where the fitted curves are plotted without the 
wide-band ~ 1// contribution. 
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expected noise sources. In particular, the fundamental hypothesis is that the noise is a 
(quasi-) stationary stochastic process with Gaussian statistics. The system statistics is a 
crucial issue for GW detectors and it can cause a dramatic decrease of the effective duty 
cycle. Indeed, it is safe to limit the analysis only to the periods when the experimental noise 
is well modeled, indicating that the detector is working properly. 

The output of PDl was recorded with an acquisition system identical to the one em- 
ployed for the ultra-cryogenic GW detector AURIGA: the data are sampled at 4.88 kHz 
and synchronized to UTC with a GPS clock. We have acquired data between May 31st 
2001 and June 20th 2001. The data acquisition was not continuous due to intentional inter- 
ruptions for diagnostic purposes and to system failures mainly originated by environmental 



temperature variations [|I9|. Manual relocking procedures require less than 15 minutes and 
the longest continuous locking period was 42 hours. The overall data recording corresponds 
to 183 hours. A few additional signals were sampled at 20 Hz for monitoring the dc signals 
from PDl and PD2, the temperature of the RC and the correction voltage fed to the PZT 
of the RC. 

The acquired data were elaborated through the same data analysis used for the ultra- 
cryogenic detector |^ . The analysis implements a Wiener-Kolmogorov (WK) filter to search 
for 5-like signals (triggers), i.e., for short bursts whose Fourier transform can be considered as 
constant over the effective bandwidth of the detector. A maximum-hold algorithm is applied 
to the data and for each trigger we estimate the time of arrival, the amplitude and the 
with respect to the expected shape. The latter discriminates a 5-like mechanical excitation 
of the bar, i.e., the GW signal, from spurious signals. The analysis also implements adaptive 
algorithms that update the parameters of the WK filter in order to follow slow drifts of the 
system. 

In order to verify the Gaussian behavior of the detection output signal, we study the 
distribution of the reduced (called xi i^i the following) of all triggers found by the data 
analysis: for a random variable with Gaussian statistics xl should follow a well-known 



distribution ||2J]. We focus on the 24 hours of data acquired during June 9th. The pumps 
that maintain the vacuum in the bar tank were switched on for one hour, between hl2.30 
UTC and hl3.30 UTC. We have vetoed the data acquired within this period in order to 
avoid the effect of the noise introduced by the pumps. We plot in Fig. ^ (top) the histogram 
of xl of all triggers with signal-to- noise ratio (SNR) greater than 4 [^, before and after 
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FIG. 4: Top: histogram of Xa of all triggers with SNR>4, before (dark gray) and after veto (light 
gray). Bottom: histogram of all triggers after veto and with 4<SNR<6. The solid line is the 
reduced distribution with 30 degrees of freedom fitted to the vetoed data with 4<SNR<6. The 
data refer to June 9th 2001, and the veto is applied when vacuum pumps are on. The binning 
corresponds to intervals of 0.005 for Xa- 

applying the veto (respectively, dark and light gray histogram). The tails at higher values of 
are very efficiently cut away by considering only those triggers with 4<SNR<6 (bottom 
histogram). The number of degrees of freedom used to compute Xa is 30. We can well fit 
the data having 4<SNR<6, after veto, with the theoretical distribution of Xa- the reduced 
of the fit is 0.90. 

In Fig. ^ we plot the Xa of triggers with SNR>4 versus the SNR. Most of them 
concentrate in the region of low SNR and Xa < 2. A small fraction (about 0.05%) of triggers 
follows a linear law in the log-log plot of Fig. ^. Indeed, it has been shown that a 
quadratic scale law of Xa versus the SNR is expected for signals which are not matched by 
the filter. 92% of triggers with SNR> 6 surviving the veto can be rejected as they have a 
Xa > 2.1, a threshold that corresponds to a confidence level of 3.9 x 10""^ for our 30 degrees 
of freedom. 
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FIG. 5: Scatter plot of Xa versus SNR for all triggers after veto and with SNR>4 (gray stars). The 
closed circles correspond to the vetoed triggers. The data refers to June 9th 2001, and the veto is 
applied when vacuum pumps are on. 



An independent test of the Gaussian character of the system is the distribution of the 
SNR. In this case the analytical formula for the expected distribution involves the calculation 
of integrals that are not easily solvable. We therefore compare the experimental distribution 
of SNR with that obtained by simulating with Monte Carlo methods a Gaussian system 
having the same parameters as ours (namely, frequencies, bandwidth and ratio between the 
height of the mechanical mode peaks and the background noise). The simulation output 
is passed through the same analysis as the real data, so that any deviation of our system 
from a Gaussian behavior would appear as a difference in the SNR distribution between 
the simulated and the true system. The results are plotted in Fig. ^ and the agreement is 
excellent, above all considering that 92% of the triggers at SNR> 6 are rejected by the 
test. The remaining signal distribution perfectly corresponds to the expected output of a 
system with Gaussian input noise, without any excess trigger. 
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FIG. 6: Histogram of the SNR of all triggers with SNR>4, before (black) and after veto (light 
gray). The data refers to June 9th 2001, and the veto is applied when vacuum pumps are on. The 
solid line is the distribution predicted by a numerical simulation of a Gaussian process with the 
same parameters as the detection system. The binning corresponds to intervals of 0.2 for SNR. 

D. Performance as GW detector 



The sensitivity of the bar detector is determined from Shh, defined as the spectral density 
of the total noise referred to the detector input and calibrated in terms of GW amplitude. 
Figure |^ shows y/Shh in a neighbourhood of the two modes '+' and '— '. As expected, the 
best sensitivity is achieved in correspondence of the modes, with a peak at the mode ' — ', 
that has a higher Q. 

A critical parameter for a GW resonant detector is the detection bandwidth, which is 
particularly significant for the temporal definition of the candidate events and for the syn- 
chronization of several detectors . The spectral sensitivity of the operating cryogenic bar 
detectors is characterized, with only one exception p^, by narrow quasi-Lorentzian peaks, 
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FIG. 7: Sensitivity of the GW detector, expressed in terms of the equivalent strain noise at the 
detector input. The data correspond to a one-hour average. 

due to the relatively large amplifier noise. In that case, the detection bandwidth is of the 
order of few Hz. For our detector, the sensitivity is relevant in the whole frequency inter- 
val between the modes, as expected for an optimized resonant transducer. The significant 
bandwidth exceeds the modes splitting. The output spectrum is asymmetric, but we remark 
that at 10 dB from the minimum the width is about 50 Hz. 

The amplitude of a GW burst that would be detected with unitary SNR is /imin = 
3 X 10~^^. This is equivalent to a standard pulse of O.SMqC^ converted into GW at the 
distance of 10 kPc, i.e., at the galactic center. This performance is not of astrophysical 
interest, as expected for the room temperature operation, but it is very promising: once the 
system is upgraded both in the optics and in the detector design to operate at ultra-cryogenic 
temperatures, the sensitivity is expected to increase up to /imin ~ 10~^°. 

IV. CONCLUSIONS 

Wc have operated a room temperature GW bar detector equipped with an opto- 
mechanical readout. The sensitivity is limited by thermal noise, due to the high temperature 
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and low mechanical quality factors. Both parameters will improve in the case of cryogenic 
operation. The statistics of the detector noise is stable and Gaussian as expected. The sensi- 
tivity is enough to evidence the failure of the normal mode expansion due to the presence of 
inhomogeneously distributed losses. The achieved results are a significant step towards the 
realization of an ultra-cryogenic resonant detector equipped with optomechanical readout. 
This is predicted to improve the bandwidth of the AURIGA detector by more than one 
order of magnitude with respect to its present status. 
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